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A novel regime is proposed where, employing linearly polarized laser pulses at intensities 
10 21 Wcm -2 as two order of magnitude lower than earlier predicted [T. Esirkepov et al., Phys. 
Rev. Lett. 92, 175003 (2004)], ions are dominantly accelerated from ultrathin foils by the radiation 
pressure, and have monoenergetic spectra. In the regime, ions accelerated from the hole-boring 
process quickly catch up with the ions accelerated by target normal sheath acceleration (TNSA), 
and they then join in a single bunch, undergoing a hybrid Light-Sail/TNSA acceleration. Under 
an appropriate coupling condition between foil thickness, laser intensity and pulse duration, laser 
radiation pressure can be dominant in this hybrid acceleration. Two-dimensional PIC simulations 
show that 1.26GeV quasimonoenergetic C 6+ beams are obtained by linearly polarized laser pulses 
at intensities of 10 21 Wcm -2 . 

PACS numbers: 52.38.Kd, 52.50.Jm, 29.25.-t, 52.65.Rr 



Laser-driven ion acceleration from solid foils has attracted great interest due to its many prospective applications 
P, [2] including isotope production, tumor therapy, ultrafast radiography, and inertial confinement fusion. Several 
of these applications require a high-energy ion beam with a large particle number and a monoenergetic spectrum. 
Currently two main acceleration mechanisms have been identified and widely investigated: target normal sheath 
acceleration (TNSA) 0, Q and radiation pressure acceleration (RPA) [El-Q. In TNSA, the acceleration of ions is 
due to the strong electrostatic sheath field at the rear of the foil, created by hot electrons generated from the foil 
front side via the oscillating j x B heating of linear polarized (LP) lasers. The ion beams produced by TNSA are 
typically characterized by low particle density, large divergence and broad energy spread. RPA in principle is a very 
efficient scheme for obtaining high-quality monoenergetic ion beams, via the so-called "hole-boring" (HB) ^ and 
"light-sail" (LS){zl, Q stages. It is generally accepted that an effective implementation of the RPA scheme at realistic 
irradiances is facilitated by the use of circularly polarized (CP) laser pulses 0, H, 0|, as this will provide a constant 
ponderomotive drive. This non-oscillating ponderomotive force applies a steady pressure to the foil front surface and 
piles up electrons into a compressed layer, inducing an intense charge separation field that accelerates the ions. This 
process, known as HB, is repeated till the laser pulse punches through the foil. At later times the compressed electron 
and ion layers combine together to form a quasineutral plasma slab undergoing LS RPA [8[. 

By contrast, for LP lasers the oscillating j x B term causes heating and recirculation of hot electrons through the 
target, resulting in a large sheath field at the rear side of the foil. The ions at the rear are accelerated by TNSA 
and the foil undergoes rapid decompression, breaking the equilibrium condition Q required for LS RPA. A RPA 
dominated regime ("laser-piston") has been identified by Esirkepov et al [lo| at the extremely high intensity of 

1 > 10 23 Wcm -2 . The key principle is using the ultra-strong radiation pressure to accelerate protons to relativistic 
energy during few laser cycles before the rapid growth of the rear sheath field and the heavy decompression of the 
foil, so that the protons can quickly catch up with the electrons in an accelerating plasma slab. Recently, Zhuo et 
al pjj , discussed a dual-peaked electrostatic field acceleration for production of quasimonoenergetic ion beams using 
LP pulses. In reality, the situation they discuss is qualitatively the same as the laser-piston regime in [10[, as they 
also use a very high laser intensity (5.5 x 10 22 Wcm -2 ) resulting in proton acceleration to relativistic energy within 

2 laser cycles. There is however a clear interest in identifying regimes where RPA can dominate the acceleration at 
intensities more reachable with current or near term laser systems employing LP laser pulses. 

In this letter, we identify a new RPA dominated regime of quasimonoenergetic ion beam generation for LP laser 
pulses at intensities of 10 20 — 10 21 Wcm -2 . In this regime, with an appropriate coupling condition between foil 
thickness, laser intensity and pulse duration satisfied, RPA can be dominant during the competition between RPA 
and TNSA mechanisms in a hybrid acceleration stage. Such a coupling condition has been analytically derived and 
verified by two-dimensional (2D) particle- in-cell (PIC) simulations. These show that 1.26GeV quasimonoenergetic 

C 6+ 

ion beam with energy spread ~ O.lGeV and angular divergence < 5° is obtained by irradiation of 80nm ultrathin 
foils with LP pulses at intensities 10 21 Wcm -2 . As expected the simulations also show that a multispecies 



foil helps to stabilize the LS RPA component in the hybrid acceleration of lighter ions, resulting in more pronounced 
spectral peaks, similar to what observed from "pure" RPA using CP lasers. 
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When a LP laser pulse irradiates a solid foil target, the electrons in the skin depth of the front surface are accelerated 
by the ponderomotive force [l2[ 

/p = "4^ V|i?i|2(1 " COS2 ^ ) ' (1) 

which includes non-oscillating and oscillating terms. In Eq. (pQ), El is the amplitude of the laser field oscillating at 
frequency cj, m e and e are electron mass and charge, t is the time. Electrons are pushed inward by the non-oscillating 
ponderomotive term and pile up in a compressed layer, inducing a charge separation electric field E z which accelerates 
ions in a pistonlike manner, in the process known as HB RPA [6, 7, 9]. From the momentum balance, the HB velocity 
is obtained as |6| 



Vb ^n c m e /njmia 
c 1 + \f n c m e /riimia ' 

where rii and mi are ion density and mass, n c = m e uj 2 / e 2 is the critical density, a = e\EL\/mcou is the laser normalized 
amplitude and c is light speed. Assuming that the laser field is perfectly reflected, the ions can obtain the maximum 
velocity Vij ront = 2v h /(l + v 2 /c 2 ). 

However, due to the oscillating ponderomotive term, a fraction of the electrons in the skin depth are heated 
rapidly by the j x B effect and travel through the target at speed close to c. A hot electron cloud forms at a 
distance of the Debye length A^ = (kBTh/^irrihe 2 ) 1 / 2 away from the target rear edge, which induces a sheath electric 
field accelerating ions according to TNSA mechanism. This sheath field scales as the hot-electron temperature Th 
multiplying the effective density in the cloud n^, i.e., E s h ea th ~ (^rih Th) 1 / 2 • Based on the expanding plasma model, 
the maximum ion velocity obtained by TNSA can be estimated as [13|, LL4[ 

Vi.rear = C h [2 ln(u pi t) + In 2 - 1] , (3) 

where Ch = (Th/rrii) 1 / 2 and uj p i = (AirnhZe 2 /mi) 1 / 2 are the hot electron sound speed and the local ion plasma 
frequency. Furthermore, for a solid foil target of steep density gradient, Th is obtained as [15] Th = (7n c /no e ) 1 ^ 2 [(l + 
a 2)i/2 _ i] meC 2 ? w here 7 = (1 + a 2 ) 1/2 . 

If the foil thickness is less than half of the pulse length, i.e., Iq < ctl/2, but thick enough to prevent laser from 
punching through, it has been shown [l6| that the recirculation of hot electrons through the target increases the 
effective hot electron density, resulting in enhanced rear sheath field and ion TNSA. The hot-electron recirculation 
also leads to rapid foil deformation, breaking the equilibrium condition required for LS RPA. Therefore, TNSA 
usually dominates for ion acceleration by LP laser pulses. The ion acceleration in the sheath field will end achieving 
the terminal velocity of Eq. (j3j) after the laser pulse is over. 

To achieve RPA dominance with LP laser pulse, a first condition is to use a foil target sufficiently thin that the 
hole-boring can reach the foil rear surface early enough within the laser pulse duration, that is, Iq < v^tl- Under this 
condition, the ions from the foil front side accelerated by HB RPA combine together with those undergoing TNSA at 
the foil rear and experience a hybrid TNSA-LS acceleration. Furthermore, to ensure that RPA dominates over TNSA 
in the hybrid acceleration, the condition Vij ront > Vi yrear should be satisfied at the time thb = lo/vb when HB reaches 
the foil rear surface. This indicates that the ions accelerated by HB RPA can overtake the TNSA- accelerated ions in 
the LS stage. From Eqs. ((2]) and ((3]), we obtain that the coupling condition of the RPA dominance is 

1 n c l 1 PnT r IZn c a n u/4 l-ln2 

a < T < aexp W ==( 7 4 

7r no A 27ry7no V ^0 v7 — 1 7 n c 2 

where no is initial foil electron density and rth = jn c is assumed to be the effective hot electron density; A is laser 
wavelength; A and Z are ion mass and charge number satisfying no = Zrii and mi = im p . The left part of the 
condition (|4]) satisfies the requirement to avoid complete electron blown-out from the foil target [H, [s[ . 

If the condition ((4]) is satisfied, after the HB ends the ion beam acceleration is dominated by RPA in LS stage 
where TNSA contributes only as a small leakage from the rear surface of the accelerated slab. The peak velocity of 
ion motion in the LS can be described [f| 0] as 

dvj = 1 \E L \ 2 1 - Vj/c 
dt 27rriimilo 7 3 l + Vi/c 

where \El\ = (m e ccj/e)a. According to Eqs. (|2])-(|5j), Fig. [T]plots three possible cases of ion hybrid acceleration from 
foils by LP laser pulses. If an ultrathin foil is used with the condition (|1J) satisfied, ion acceleration is in the RPA 
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FIG. 1: Three regimes for ion acceleration from thin foils by LP laser pulses, based on Eqs. J2])-®: RPA- dominated (blue 
and purple lines), competing of RPA andTNSA (black and purple), and TNSA-dominated (purple and green), where the cross 
point corresponds to the transition time from HB RPA to LS RPA. 



dominated regime (blue line), where high-energy quasimonoenergetic ion beams can be obtained. If the foil is thicker 
with the condition ^ violated but thin enough to satisfy Iq < v tl, ion acceleration takes place in a hybrid stage 
where RPA and TNSA compete (black and purple lines). Since the energy increase of RPA grows more rapidly than 
TNSA, quasimonoenergetic ion beams may be obtained after a long acceleration time, that is, a long pulse duration tl 
is required. If the foil is further thicker with Iq > v tl-, ions undergo pure TNSA acceleration (purple line). Therefore, 
quasimonoenergetic ion beams can be obtained by using either an ultrathin foil target or a long laser pulse even at 
intensities much lower than those discussed in lol Hl|. 



In order to verify the different regimes discussed above, 2D PIC simulations are run with the code "ILLUMINA- 
TION". We choose a LP laser pulse with I = 3 x 10 21 Wcm- 2 (a = 47), A = 1.0/xm and super-Gaussian intensity 
distribution exp[— (r/ro) 4 ] propagating along z-axis, where ro = 5/im is the spot radius. The laser pulse has a trape- 
zoidal temporal profile of duration 20to (to = 27t/cj), consisting of a plateau of 18to and rising and falling times of lro 
each. The fully-ionized carbon foil with electron density no = 600n c , C 6+ density ni = 100n c and thickness /o = 80nm 
is located at z = 3/im, where the condition (j4]) is satisfied. In the simulations, 5000 cells along z-axis and 3072 cells 
transversely along x-axis constitute a 20 x 24.576/im simulation box. Each foil cell is filled up with 576 electrons and 
96 C 6+ ions. 

Figures Efa)-|2fd) show density maps of electrons (n e ) and C 6+ (n^) at time t = 15, 23, 32 and 53fs. Part (a) 
corresponds to the initial HB stage. On the one hand, due to the j x B heating, the hot electrons of temperature 
Th = (7n c /noe) 1 ^ 2 [(l + a 2 ) 1 / 2 — l]m e c 2 « 7MeV quickly travel through the target and form a hot electron sheet 
of density about lln c at a distance of their Debye length ~ 0.2/im away from the foil rear surface. From Fig. 
[2fe) we see that a strong sheath electric field E s heath — 6m e cuj/e = 1.93 x 10 13 V/m is induced which accelerates 
ions from the foil rear side. On the other hand, electrons at the front side are pushed inward by the non-oscillating 
radiation pressure and pile up in a compressed layer, leading a much more intense charge separation electric field 
E separation = 22m e cuj/e = 7.0 x 10 13 V/m, seen in Fig. Efe), which accelerates ions from the foil front side. 

Figure Efb) corresponds to the transition time when HB reaches the foil rear surface. We see that the compressed 
electron layer catches up with the rear hot electron sheet; correspondingly the ions that undergo HB RPA also keep 
up with those at the rear which experience TNSA because E separat i on > E s heath- At this stage, a majority of ions 
are bunched together undergoing RPA, but they are preceded by a small number of TNSA ions. This ion bunch is 
accelerated by a combination of the charge separation field at its front and the sheath field at its rear, see Fig. [2ff), 
undergoing a hybrid LS-TNSA acceleration. Since the condition (j4]) is satisfied, E separation remains always larger 
than E s heath, resulting in quasimonoenergetic ion beams in Fig. Efc). This hybrid acceleration dynamics can also be 
seen in the ion phase plots [2fg) and[2fh): a majority of ions undergo typical RPA "rolling- up" and concentration, 
whereas a small number of ions precede like a leakage due to TNSA. 

Figures [3f a) and[3fb) show both a plot of C 6+ density at £ = 83fs when the laser pulse ends and its energy spectra at 




3 3.2 3.4 3 , x 3.5 3 3.4 3.8 3 3.5 4 4.5 5 

z(jim) z(jxm) z (jim) z (fim) 




FIG. 2: (a)-(d) Electron (n e ) and ion (C 6+ , m) densities at respectively £=15, 23, 32 and 53fs for a 80nm fully-inoized carbon 
foil irradiated by LP laser pulses at intensity Io = 3x 10 21 W/cm 2 , wavelength A = 1.0/xm and duration tl — 20ro = 67fs, where 
the foil electron density is 600n c ; (e) and (f) the longitudinal profile of m and electrostatic field E z at t = 15 and 23fs, which 
shows ions are accelerated by a hybrid effect of the charge separation field and the rear sheath field, i.e., undergoing hybrid 
RPA-TNSA; (g) and (h) C 6+ phase space distributions at t = 23 and 32fs for particles with \x\ < 2/rni, which show the typical 
characteristics of hybrid RPA-TNSA mechanism. 



£ = 53, 68 and 83fs. We see that a quasimonoenergetic C + beam with a well defined peak energy of 1.26GeV, density 
4n c and divergence < 5° is obtained at t = 83fs by the dominant RPA. The low-energy trailing part in the energy 
spectrum is originated from an extended area (wings) of the foil around the laser focal spot, where the laser field 
is weaker with the condition (4) violated and the transverse instabilities [ItJ grows up quickly, so most likely TNSA 
dominates. Figs. [3fc) and [3fd) show plots for a foil thickness of 200nm where condition (j4]) is violated. We can see 
that the energy spectrum shows no pronounced peak because RPA is not dominant in the hybrid acceleration. 

In order to emphasize the practical utility of this new regime, we present the scaling of the required foil thickness 
lo with the laser power in Fig. EJe), where the laser pulse duration of 66fs and wavelength 1.0/xm are fixed, and 
the same foil target density as above is chosen. We can see that the RPA dominated regime can be achieved by 
using ultrathin foils (10s of nm) on the contemporary high-contrast 100s of Terawatt and/or near-term Petawatt 
laser systems. Our recent experiment [2l[ has also observed narrow band features in both proton and C 6+ spectra 
from lOOnm copper foil by Vulcan laser (energy 200J, duration ~ 500fs) at Rutherford Appleton Laboratory, which 
has further verified the regime. Furthermore, note that the trapezoidal temporal profile of laser pulse chosen in the 
simulations are only intended to to show the exact validity of the condition (4) as the latter is derived by assuming 
a is a constant. The regime itself in principle is not sensitive to laser pulse profile, which only requires the general 
condition Vij ront > Vi^ rear [from Eqs. (2) and (3)] to be satisfied. 

A typical experimental target is always a solid foil contaminated with lighter ions (carbon and hydrogen). As 
expected [sKl! 20], a multispecies target helps to stabilize LS RPA of the light ion species so that a more pronounced 
peak in the energy spectrum can be achieved. Fig. |4]shows density map and energy spectrum of the acceleration case by 
choosing a fully-ionized copper foil with carbon contaminants, where the density ratio between Cu 24+ and C 6+ species 
is 5 : 1. Their total mass densities and the other parameters are the same as above. We see that a more pronounced 
monoenergetic C 6+ beam is obtained with density 3n c , divergence < 5° and higher peak energy 1.44GeV (120MeV/u), 
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FIG. 3: (a) C 6+ density at t — 83fs when the laser pulse ends for the acceleration of Fig. [2j (b) C 6+ energy spectra at t = 53 
(black), 68 (blue) and 83fs (red) for Fig. [2j which have a well-defined peak as RPA is dominant in the hybrid acceleration, (c) 
and (d) are those for the case by increasing foil thickness to /o z =200nm where condition Q is violated, showing no pronounced 
peak as RPA is not dominant, (e) shows the scaling of the required foil thickness for the RPA dominated regime with the 
laser power by fixing laser pulse duration at 66fs and wavelength 1.0/im, where the target density are the same Fig. [2] The 
corresponding points for the spectral cases of (b) and (d) are also respectively marked in (e). 



whereas the Cu 24+ species has insufficient co-moving electrons and undergoes Coulomb explosion resulting in a 
broadened energy spectrum. Note that here the particle number (dN/dE) is about 10 10 MeV _1 , much larger than 
peaks obtained from the improved TNSA 

To summarize, we have reported a new RPA dominated regime of quasimonoenergetic ion beam generation from 
ultrathin foils for LP laser pulses at intensities two-orders of magnitude than earlier predicted [10| . A criterion for this 
regime has been identified analytically and verified by 2D PIC simulations, which indicates that the RPA dominated 
regime can be achieved with LP laser pulses by using either a sufficiently thin foil target or a long enough laser pulse. 
It is shown by 2D PIC simulations that quasimonoenergetic C 6+ ion beams with peak energy about 1.26GeV, particle 
number 10 11 and beam divergence < 5° are obtained from 80nm carbon foils by LP laser pulses at intensities of 
10 21 Wcm~ 2 . 

We acknowledge helpful discussions with A. Macchi, J. Schreiber and B. Dromey. 
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FIG. 4: Ion acceleration from a 80nm fully-ionized copper-carbon mixed foil target by LP laser pulse, where the density ratio 
between Cu 24+ and C 6+ is 5 : 1. The total mass density and other parameters are the same as those in Fig. [2] (a) is the plot 
of both Cu 24+ and C 6+ densities at £ = 83fs when the laser pulse ends; (b) shows the corresponding energy spectra. The dashed 
blue and green lines are the energy spectra for C 6+ in the cases of thicker foil with Iq — 200 and 800nm. 
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